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Abstract 

We use the fermion zero-modes in the background of multi-caloron solutions with non- 
trivial holonomy as a probe for constituent monopoles. We find in general indication 
for an extended structure. However, for well separated constituents these become 
point-like. We analyse this in detail for the SU(2) charge 2 case, where one is able 
to solve the relevant Nahm equation exactly, beyond the piecewize constant solutions 
studied previously. Remarkably the zero-mode density can be expressed in the high 
temperature limit as a function of the conserved quantities that classify the solutions 
of the Nahm equation. 
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1 Introduction 



To describe regular monopoles in gauge theories a Higgs field is required. This de- 
fines the abelian subgroup of the gauge field. Yet in the full non-abelian theory there 
is no Dirac string and a regular solution results, the well-known 't Hooft-Polyakov 
monopole 0. In the strong interactions no such Higgs field should be present, but 
nevertheless it has been conjectured that a dual superconductor description, in which 
monopoles form the dual charges that condense, could explain confinement [2]. This 
scenario receives some support from the studies in supersymmetric theories through 
Seiberg-Witten duality [3 , although also the old center vortex picture is still under 
active study [I]. Lattice studies based on abelian and center projections, and their 
respective notions of dominance |S] are the main means through which one tries to 
address these issues. One relies on so-called gauge fixing singularities to identify the 
relevant monopole [S] or center vortex degrees of freedom. 

A more recent alternative to study the monopole content of gauge theories, without 
the need of addressing singular configurations, gauge fixing, nor introducing an extra 
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Higgs field, has been through calorons, which are instantons at finite temperature. It 
has been found that calorons are actually made up from constituent monopoles [11 El El , 
which becomes most apparent when the background Polyakov loop is non-trivial (as in 
the confined phase), and the size of the caloron is larger than the inverse temperature 
(the extent of the euclidean time direction). The background Polyakov loop is defined 
in the periodic gauge A^{x, t) = A^(x,t + f3) by its asymptotic value, or the so-called 
holonomy, 

rP 

Voo = lim Pexp( / A (x, t)dt) = # T exp(27udiag(/ii, /i 2 , • • • , Vn))g, (1) 

x >oo Jq 

where g is the gauge rotation used to diagonalize Voo, whose eigenvalues exp(27ri/Xj) 
can be ordered on the circle such that /ii < /i2 < . . . < fi n < AWi, with fj, n+ j = 1 + fij 
and Yn=i Hi = 0. The masses of the constituent monopoles are given by 87r 2 z/j//3, with 
Uj = fij+i — fij, which add up to 8n 2 //3, consistent with the instanton action. 

Solutions are known explicitly [TD] for SU(n) charge 1 calorons. The n constituents 
can have any spatial location, although all have the same location in time (they do, 
however, become static when well separated), and n— 1 abelian phases complete its An 
parameters. Charge k calorons can be viewed as composed of kn monopoles, of which 
a class of axially symmetric configurations was constructed explicitly [TTj . 

The purpose of this paper is to study in more detail these higher charge calorons, 
where the emphasis is on constructing the chiral fermion zero-modes. Charge 1 calorons 
have exactly one fermion zero-mode, which was shown for well separated constituents 
to be supported on one and only one constituent ^21 [HI! • We may change the con- 
stituent that supports the zero-mode, by changing the fermion boundary conditions 
from (anti)periodic, to being periodic up to a phase exp(2iriz) (from now on we will 
use the classical scale invariance to set (3 = 1). For z G the zero-mode is 

localized to what we will call type j constituent monopoles (with mass 87r 2 ^-, and the 
appropriate f/ n_1 (l) charge associated to their embedding in SU(n)). 

Lattice evidence has been gathered over recent years that these monopole con- 
stituents are present in dynamical configurations in the confined phase of gauge theo- 
ries for SU{2) using cooling [HI EH], an d f° r SU(3) using fermion zero-modes [THj as a 
filter. It is somewhat of a puzzle that these constituent monopoles had not been seen 
in earlier cooling studies (apart from when using twisted boundary conditions [T7j). 
That they remained unnoticed when using fermion zero-modes as a filter is, however, 
simply a consequence of the fact that these studies were restricted to the use of fixed 
fermion boundary conditions. Only when cycling through boundary conditions spec- 
ified by periodicity up to a phase, the SU(3) charge 1 instanton configurations will 
show three separate constituent monopoles. In Fig. ^ we show a typical example based 
on the exact solutions for SU(3), closely following the observed behavior [TE] based 
on actual lattice simulations in the confined phase, which guarantees the background 
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Polyakov loop to be non-trivial. In the high temperature phase, where the Polyakov 
loop is trivial, two of the constituents are massless and only one peak will be seen. 
These massless constituents are interesting in their own right, giving rise to so-called 
non-abelian clouds [IB], but they will not concern us here. 




Figure 1: The logarithm of the properly normalized zero- mode density for a typical 
SU(3) caloron of charge 1, cycling through z. Shown are z = fij (for linear plots see 
Fig. EJ) and three values of z roughly in the middle of each interval z G \p,j, fij+i]- All 
plots are on the same scale, cutoff for values of the logarithm below -5. The zero-mode 
with anti-periodic boundary conditions is found at z — 30/60. For the action density 
of the associated gauge field, see Ref. [T3"| IT$|. 




Figure 2: The properly normalized delocalized zero-mode densities for z = fij, on a 
linear scale (cmp. Fig. QJ. 

We do not only study the fermion zero-modes for the higher charge calorons to 
compare with lattice simulations, but also as a tool to understand to which extend the 
constituent monopoles can be unambiguously identified in the caloron solutions. In the 
high temperature limit the non-abelian cores of the monopoles shrink to zero-sizes, and 
one is left with abelian gauge fields. Without taking the high temperature limit, but 
excluding the non-abelian cores, the same physics describes what we called the far field 
region [H]. It would be desirable if the abelian field in this region is described by point- 
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like Dirac monopoles (actually dyons because of self-duality), when the constituents 
are well separated. For charge 1 calorons and the class of axially symmetric solutions 
studied before ^l] the density of the fermion zero-modes become Dirac delta functions 
at the locations of the constituent monopoles in the high temperature limit, for any 
constituent separation. This infinite localization in the high temperature limit can 
be understood from the fact that for most z values there is an effective mass for the 
fermions. Therefore, by studying these zero-modes in the general case, we hope to 
learn more about the localization of the monopoles. 

It is not directly obvious that any higher charge caloron can be described by point- 
like constituents in the far field limit. The tool to construct these solutions is the Nahm 
equation [201 Ej, which is a duality transformation that maps the problem of finding 
charge k calorons to that of U (k) gauge fields on a circle with specific singularities at 
z = fij. In the cases studied so far, the dual (Nahm) gauge field can be made piecewize 
constant, from which one easily reads off the constituent monopole locations. In gen- 
eral, however, the Nahm gauge fields depend non-trivially on z, and it is important to 
understand what this implies for the localization of the constituents. Quite remarkably, 
we will nevertheless find that in the high temperature limit the fermion zero-mode den- 
sity does not depend on z and can be expressed in terms of the conserved quantities 
of the Nahm equation. We compute it explicitly for charge 2, revealing in general an 
extended structure. However, the extended structure collapses to isolated points for 
well separated constituents. This is so as long as z ^ fij, which is the value where the 
localization of the zero-modes jumps. We separately study in the high temperature 
limit the case where z = fij, for which the fermion zero- modes delocalize (decaying 
algebraically). Support of the zero- modes is now on constituent monopoles of types 
j — 1 and j. We analyse this in detail for the class of axially symmetric solutions of 
Ref. JXTJ, in the light of some puzzles concerning so-called bipole zero-modes [2~Tj . 

This paper is organized as follows: first we describe the caloron zero-modes in 
Sect. E] and set up the Green's function calculation in Sect. H2 to be able to discuss 
the zero-mode and far field limits. Sect. H] deals with the case where z = fij, for which 
zero-modes are delocalized. In Sect. El we relate the zero-mode density for z ^ fij to 
conserved quantities of the Nahm equation and study in Sect. [HI the general solution 
for SU (2) charge 2 caloron. We conclude with a discussion of the implications of our 
results and the problems that still need to be addressed. Two appendices provide the 
details for the zero-mode and far field limit calculations. 

2 Fermion zero-modes 

We wish to construct the zero-modes of the chiral Dirac equation in the background of 
a self-dual gauge field at finite temperature. The Dirac equation in its two-component 
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Weyl form, with oy, = a t = (1 2 , — zt) (Tj are the usual Pauli matrices), reads 

Di> z (x) = a^D^^x) = + A^x))^ z {x) . (2) 

Assuming the gauge field is periodic in the imaginary time direction, with period (3 = 1, 
we seek the zero-modes that satisfy the boundary condition 

i> z {t + l,x) = exp(-2mz)if z (t, x), (3) 

A simple abelian gauge transformation, ^ z {x) = exp(2iTizt)^/ z (x) , makes the zero- 
mode periodic. This gauge transformation replaces the gauge field by A^(x) — 2niz tl t n , 
but it does not change the field strength, such that existence of the appropriate number 
of zero-modes is guaranteed by the index theorem. 

The caloron solutions are obtained by Fourier transformation, reformulating the al- 
gebraic ADHM (Atiyah-Drinfeld-Hitchin-Manin) construction [22] of multi-instanton 
solutions in R 4 , as the Nahm transformation For this the instantons in R 4 period- 
ically repeat themselves in the imaginary time direction up to a gauge rotation with 
Voq. This Fourier transformation also selects out of the infinite number of fermion 
zero-modes in R 4 , those that satisfy the correct periodicity. 

In the ADHM formulation the k normalized fermion zero-modes for a charge k 
SU (N) instanton are given by (e = IT2 = <7 2 , i = 1, . . . , n is the color index, I = 1, . . . , k 
the charge index and J = 1 , 2 the spinor index) 

¥ u (x) = vr 1 (<p- l / 2 {x)u\x)f x e)\ i , <f>{x) = t n + u\x)u(x), (4) 

with v)(x) given explicitly in terms of the ADHM parameters by 

u\x) = X(B - l k x)'\ B = a^B^ x = x^a^, (5) 

where A = (Ai, ...,A^), with A| a two-component spinor in the n representation of 
SU(n) (A can be seen as a n x 2k complex matrix), and B^ hermitian k x k matrices. 
As 4>(x) is an n x n positive hermitian matrix (for n = 2 proportional to I2), its square 
root is well-defined. We also recall the gauge field is given by 

A^x) = <j)- 1/2 (x)(^(x)d^u(x))(f ) - 1/2 (x) + 1/2 (x)^0- 1/2 (x). (6) 

which can be shown to be self-dual provided the quadratic ADHM constraint is satisfied, 

tf\ + {B-± k x)\B-± k x)=a f-\ (7) 

which implicitly defines f x as a hermitian k x k Green's function, thereby completing 
the description of Eq. (@J). A further simplication 0122] will be helpful, which uses the 
fact that 2{B^ — x^) = d^f' 1 and vf(x) = <p(x)\f x (B — tkx)\ implying 

A,(x) = ^ 2 (x)Xf}^dJ x XU 1/ \x) + ^~ 1/ \^):d^ 1/ \x)l 
<F\x) = 1-A/.A*, ¥ iT (x) = (2n)- 1 (^ 2 (x)Xd fl f x a lx e) l iI , (8) 
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with r]^ = rf^Gj = o\ il (Jv\ the anti-selfdual 't Hooft tensor. 

As mentioned above, calorons are obtained by arranging the instantons in R to 
be periodic (up to a gauge rotation). The time interval [0,1] will contain as many 
instantons as the topological charge of the caloron. One splits the charge index I as 
I = pk + a, where a labels the k instantons in the interval [0, 1] and p labels the infinite 
number of repeated time-intervals, playing the role of the Fourier index. We find, 
suppressing the gauge and spinor index (cmp. Ref. [T2"t H3]). 

9%x) = (27r)-V 1/a (aO0„ f dz'X b (z')t(z', z)a,e. (9) 

J 

where the Fourier transforms of A and f x are denoted by X a (z) and f ba (z',z). The 
fermion zero-modes thus constructed are in the so-called algebraic gauge, for which 
iff z (t + 1,2?) = exp(— 27riz)V O0 J & z (t, x). In this gauge all components of vanish 
at spatial infinity and the non-trivial holonomy is encoded in the boundary condition, 
which for the gauge field reads A^(t + 1, x) = VooA^it, x)V^ . A simple time-dependent 
gauge transformation allows one to transform to the periodic gauge, after which Aq 
goes to a constant at spatial infinity. 

To encode the appropriate periodicity in the ADHM parameters we need to take 
Apfc+a = 'P^oCa [TTj . Introducing the projections P m on the eigenvalues of Poo, i.e. 
Voo = ELi exp(27rz/i m )P m , we find X a (z) = £^ =1 5(2 - yU m )-P m Ca, which makes the 
expression for the zero-modes particularly simple, 

n 

& z (x) = (2vr)- V 1/2 (x) £ P m Cb^edJ^(fx m , z). (10) 

m=l 

The fact that we are dealing with higher charge calorons, is reflected in the presence 
of the indices a,b = 1, . . . , k. Making use of the well-known identity [23 I2H 122] 
in i? 4 , ^ l iI (x)*^^(x) = —{2T[)^ 2 d'j L f l J n ) Fourier transformation gives the appropriate 
expression for the caloron zero-mode density 

¥ z {x)^{x) = -{2^)-Xt{z,z). (11) 

Using the fact that \\m\^ 00 \x\f^ b {z, z) = irS ab , the zero-modes ^ a z are seen to be 
orthonormal. 

We close this section by remarking that for 577(2) an alternative construction is 
possible, as part of the Sp{n) series (since Sp(l) = SU(2)). The ADHM construction 
for Sp(n) is based on quaternions. In particular A; is assumed to be a quaternion, 
whereas 7? M is now real-symmetric. All formulae presented above remain valid, but it 
should be noted that the transformation A — > XT\ — > TB^T^, with T £ U(k), 
leaving the gauge field and the ADHM constraint untouched, has to be replaced by 
TeO(k). 
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3 Zero-mode and far field limit 



As we have seen above, all physical quantities can be reconstructed, once we have 
found the Green's function f® \z,z'). Here we review the necessary ingredients. We 
start with the fact that the Green's function is defined through an impurity scattering 
problem ^1] 

{-^ + V(z;x)\f x (z,z')=47r 2 ± k 5(z-z'), (12) 

where f x (z, z') is related to f x (z, z') by a U{k) gauge transformation 

f x (z, z) = g(z)f x (z, z')g\z ), g(z) = exp (2m(£ - x l k )z) . (13) 

The "potential" V, which includes "impurity" contributions, is determined by the 
(dual) U(k) Nahm gauge field A^(z) 

V(z;x) = 4n 2 R 2 (z;x) + 2vr^5(z - p m )S m) S m = g{p m )S m 9 ] {p m ), 

in 

R 3 (z; x) = Xj t k - {2mY x g{z)A j {z)g\z), S% = ntr 2 (CP m (b) ■ (14) 

Fourier transformation of gives (2ixi)~ 1 5(z — z'){5 lJL Q'l k -j^ + A^(z)) and defines the 
Nahm gauge field. We have further used the fact that one can choose a U (k) gauge in 
which Aq{z) = 2tii^q is constant, which itself can be transformed to zero by g(z). Note 
that g(l) plays the role of the holonomy associated to the dual Nahm gauge field A^(z). 
Crucial is that we can formulate the zero-mode and far field limits without specifying 
the solutions of the Nahm equations. These Nahm equations, which are equivalent to 
the ADHM constraint, are given by 

■j-M*) + [Mz)>M*)\ + \e M [A k (z),A,(z)} = 2m S(z - fi m )pL 
ctz m 

= -7rtr 2 (ClPmCbr) ■ (15) 



P 



ab 
m 



This will be discussed in more detail in Sect. El 

To solve the second order equation for f x (z,z') it is convenient to convert it to a 
first order equation, involving 2k x 2k matrices, 

i ( u\\( U(z, ^) \ _ _^ S ( Z - f V (16) 



\dz \V(z;x) J J \fj x (z,z')) K '\U 

which can be solved as 



£f(' Z J ) = ~^ 2 W(z, z ) {(t 2k - - 6(z' - z)l 2k } W~\z\ 

dzJ x \ ' / ' 



(17) 
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where zq can be arbitrary and 



22 / l k . 
V(z;x) 1 



21 



, ^=^(1)^ + 1,20) 



W/(z 2 , Zi) = Pexp 
In particular, one can show that JJJ 

- |Tr n F^(x) = -JflJ^ log^(x), iP = det (ie-™°(l 2 , - F zo )/V2) . (19) 

To isolate the exponential contributions in Eq. (|18jl . one introduces two solutions of 
the Riccati equation [TT] . 

it(z) a ±i^i^(z) = ^;5). (20) 
Since i?(z; x) — > xl^ for |af| — > oo, we find in this limit that R^(z) — > \x\tk- Defining 

fm( z )= Pex P ±2VT /" R^(z)dz , Z G [Hm, flm+l], (21) 

J fj,rn 

we see that /^(z) — > exp {±2 r w\x\{z — /x m )±&). These are the exponentially rising and 
decaying solutions of Eq. (fT2J) . in terms of which we can rewrite for z, z' e (/i m , fi m +i) 

W(z, z') = W m (z)W-\z'), W m (z) = W m (z)F m (z), (22) 

where F m (z) contains all exponential factors, 

*"M = U£ W -2^))' ^-(^ /-(,))• <*» 

As an illustration let us assume R(z;x) 2 = (xtk — eY m ) 2 = R 2 W independent of z, 
for z G [p, m , fJim+i], as is the case for charge one |Hj and a class of axially symmetric 
solutions [TT] (see Sect.EJ). We then find R^z) — Rm( z ) = Rm an d fm( z ) = exp(27r(,s— 
fi)R m ). Diagonalizing R m defines k locations (the eigenvalues of Y m ), which are sharply 
defined and give rise to point-like constituents in the high temperature limit. When 
R(z; x) is not piecewize constant, these locations are a priori not sharply defined. We 
still expect the cores to be the regions in x where R(z; x) remains small (cmp. Eq. ([14))) . 
The separations between the cores of monopoles of different type is controlled by the 
discontinuity in Aj(z), Eq. ()15|) . and can in general be chosen large. This allows us to 
define the zero-mode limit, where x is assumed to be far removed from any constituent 
core not of type m. In technical terms that means that R^,(z) is large for all m' ^ m. 
Accordingly, f~>(z) and f^'i 2 )' 1 are exponentially small for these values of x (cmp. 
Eq. ()21j)). Results that are valid up to these exponential correction are denoted by a 
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subscript "zm" for the zero-mode limit and "ff" for the far field limit. In the latter 
case, x is assumed to be far removed from all constituents. 

In appendix A we show that for // m < z' < z < fi m+ i 

fx i Z ) Z ) = 7r ( e m(^) — e m{ Z )^m+l) ( e m — ^m e m^m+l) {&m( Z ) — ^m^m( Z ))^m i Z )• 

(24) 

(for z' > z one uses f x (z', z) = f£(z, z')) with 

e m( Z ) = fm( Z ) fmiftm+l) > e m = ^mit^m)-, &m( z ) = fm( Z ) ) 

Z> m = Ifc — 2S~ 1 i? m _i(/i m ), = l fc — 2E~ 1 i? m (/z m ), (25) 

i? m (z) = \{RZSz) + S m EE R~(Hm) + Rt-l(Vm) + S m . 

This result is valid up to exponential corrections as long as x is far removed from all 
constituents of type m' ^ m. Note that in this limit Zf^ = l fc , however, only up to 
algebraic corrections, which is why we have kept them. In Fig. El we give the zero- 
mode densities for a charge 2 axially symmetric solution in SU(2), with well-separated 
constituents. We found no differences to 1 part in 10 6 between the exact result and 
the result obtained with the zero-mode limit, Eq. (|24p . The choice of basis for the 
zero-modes involves the gauge rotation that for each interval z G [/i m , /i m +i] identifies 
the constituent locations, cmp. Sect.|3]and Ref. |TT]. This ensures that each zero-mode 
is localized to only one constituent monopole. 




Figure 3: Zero-mode densities for a typical charge 2, £77(2) axially symmetric solution. 
For comparison the action density (cmp. Fig. 2 of Ref. JT]) is shown in the middle. 
All are on a logarithmic scale, cutoff below e -3 . On the left is shown the two periodic 
zero-modes (z = 0) and on the right the two anti-periodic zero- modes (z = 1/2). 

It is now almost trivial to read off from this the far field limit for f x (z, z), where x 
is assumed to be far from all constituents. As long as z ^ /i m and z ^ /i m +i> e m( z ) an d 
e~(z) are exponentially small, and f^(z,z) = Tre^ n (z)(e^ n )~ 1 e^ n (z)R^ 1 L (z). But using 
the definitions of e^(z) and e^(z) all exponential factors precisely cancel and in the 
far field limit we are left with 

fx( z , z ) = « R m( z ), z G (// m ,// m+ i), (26) 
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which will play a very important role in Sects. ElandlHl We can also read off the far field 
limit for f x (z,z') evaluated at the impurities // m and noting that by definition 

e m(Pm+l) = ^mi^m) = H, 

fx (/^m+i) fJ"m+i) — ^(Ife — Z m+1 )Rj;{ii m+ i) = 27r£ m 1 fl , 

/^Owm) = 7r(l fe -Z+)i?- 1 (/i m ) = 27rS m 1 , (27) 

and fxifJ'm+ufJ'm) = (as well as /f(// m) /Wi) = °> usin S fx(z,z') = fl(z',z)), 
verifying the results of Ref. [TT| . Although f x (z,z) is continuous at the impurities, 
comparing Eq. (|2fij) with Eq. (|2*7j) we see, as anticipated, that in the high temperature 
limit f x (z,z) is discontinuous at the impurities. At finite temperature, the transition 
across the impurity has a "width" inversely proportional to the temperature. 

For k = 1, where e±(z) = exp(±27r(/i m+ i - z)r m ) and e^(z) = exp(±2ir(z - fi m )r m ) 
one finds for z', z G (/i m , Hm+i) 



27rsinh (27rr m (/i m+1 - 2;' + 7 m+1 )j sinh (2vrr m (z - /x m + 7+)) 

fl m (z\z) = * ; -, L ^ , (2S 

sinh (27rr m (z/ m + 7+ + 7m+i) 



which agrees with the result of Ref. [12], where only the limit with 7+ = — \ logZ+ and 
7 m+1 = — ilogZ~ +1 neglected and z' = z was considered. We stress again that the 
presence of 7 implies a subtle algebraically decaying influence due to the constituents of 
type m — 1 and m + 1 (the influence of all other constituents is decaying exponentially, 
although this is only relevant for SU(n > 3)). In terms of constituent radii r m = \x—y m \ 
and locations y m one has 

r?— IZ/m+l Um\ + r m+l r m rv-\- \Um— 1 Vml + r m— 1 ?"m /on\ 

Z m+1 - "P -» I , J Z m - "p _ -» I ~ • (+ y J 

Wm+l Urn] \ 'm+1 t '"m |2/m— 1 Vm\ i fm— 1 ~r ?"m 

For SC/(2), with Z/m+i — Vm-ij therefore — Z m+l , or 7+ — 7 m+1 , and the influence 
of the other constituent is only felt by a renormalization of the mass z/ m , 

qtt(o\ fzmt ^ cosh (2nr m (v m + 2 7m )) - cosh {2%r m {2z - /i m - /Wi)) , , 
SU(2): f. („,,)=» rro8illh (2 Ir r ro (, m + 2 7m )) ' (3 °» 

Note that for SU{2) /xi + /i 2 = and /i 2 + /i3 = 1- This leads to the well-know result for 
the monopole zero-mode density f x m (z,z) = vrr" 1 tanh (vrr m (z/ m + 27 m )), with 2 = 
for m = 1 (periodic zero-mode) and z = 1/2 for m = 2 (anti-periodic zero- mode). 



4 Bipole zero-modes 

In this section we discuss the zero-modes in the high temperature limit for z = fi m , 
which means that Pooe' 2 ™ has one of its eigenvalues equal to 1, which leads to one of 
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the components of the fermion to become massless. Indeed, using Eqs. (|1 1113127(1 we 
find that in the far field limit 

z = p m : K(x)^ b z (x) = -(2tt)- 1 ^ , t m = ^JE^J, (31) 

decays algebraically and has support on the constituents of type m— 1 and m, as is easy 
to see for k = 1, where S m = \x — y m -i\ + \x — y m \ + \y m -i — Vml- Here we will restrict 
ourselves to SU(2), particularly interesting for the axially symmetric caloron solutions, 
since in the high temperature limit its gauge field has the form of the so-called bipole 
ansatz (see Eq. (|35(l ). which always has an integrable chiral fermion zero- mode [21]. In 
the bipole ansatz all Dirac strings have to run in the same direction, but other than 
that, the locations of the self-dual Dirac monopoles can be arbitrary. However, for the 
axially symmetric caloron solutions the constituents have to alternate between opposite 
charges on a line jllj . In this case, with the solution coming from a regular caloron, 
there are always as many zero-modes as the number of constituents with a given charge 
(equal to the topological charge k). By considering the case of solutions with topological 
charge 2, we will find the expressions for the bipole zero-mode and the extra zero-mode, 
in terms of the constituent locations only (which should be possible, since the abelian 
gauge field has this property in the high temperature limit). Remarkably, we will find 
that rearranging the order of the constituents, so as to violate the constraint coming 
from the axially symmetric caloron, the second zero-mode is no longer integrable (while 
the gauge field and bipole zero-mode remain well defined). 

A particular class of axially symmetric caloron solutions is obtained by taking 

Y m = g(z) A ^lf 9K Z ) t° be piecewize constant. This can be shown to satisfy the Nahm 
or ADHM equation when we take [TT] 

n 

Ca = Pa exp(27ua a )C, a a =^2a^P m , Tr„a a = 0, Y m = Y m e, (32) 

m=l 

where p a are positive, not to be confused with 

pt = pVexp(2«(«™ - a m ))Ay m , Ay m = Ay m e = -vrtr 2 (C t P m Cr). (33) 

For SU (2) one has Ayi = —A1/2 (for SU(n > 2) a constraint on ( is required, to 
guarantee that all Ay m are parallel). We will take e = Ay 2 /|Ay 2 | = (0,0,1) (hence 



Ay2 > 0) and Voo = exp(2mp2T3) (therefore C = \j ^2/^2)- This can always be 
arranged to be the case by a global gauge rotation, and a spatial rotation. We define 
y m , up to an irrelevant overall shift, through Ay m = y m — y m -i- This fixes Y m to be 

n 

Ym = + P 2 a yrn)5ab + l(l ~ ^ab)papb £ A V] 

3=1 

(34) 
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The constituent locations are found from the eigenvalues of Y m . Although constant, 
it is not true in general that the Y m can be diagonalized simultaneously, making this 
a non-abelian solution of the Nahm equation. Yet, as we remarked before, one easily 
computes the Green's function, since i?^(z) = R m = J (xtk — eY m ) ■ (xlk — eY m ) is 

constant in z. Using that for the axially symmetric solutions p^f = S^y m /\Ay m \, one 
shows jTT] that in the high temperature limit the (abelian) gauge field can be written 
in the form of the bipole ansatz [21] 

A,(x) = -irs^Alogtfz), 4>{x) = <M*) ee ^(flit^-t) ' (35) 

In the bipole ansatz, one splits A^(x) in an isospin up and isospin down component 
(with inverted abelian charges). However, all that concerns us here is the fact that 
for any <p, A^x) as given above is self-dual (and hence a solution of the Maxwell 
equations) provided log <fi is harmonic away from Dirac string singularities (defined by 
(f)^ 1 = 0). This always gives rise to at least one normalizable zero-mode of the chiral 
Dirac equation 

i> mI {x) = (27rp)-V~ 1/2 (^)(Ti^£)m/^log0^). (36) 

were p is simply a normalization factor. Here m corresponds to the isospin component 
that survives for z = p m (with the other component related to f x {pi, P2) vanishing in 
the far field limit, see Eq. Q). 

We now work out the explicit form of all k zero-modes for the axially symmetric 
caloron solution, showing how the zero-mode in Eq. (|36|) is recovered from these. Using 
Eqs. f)l()p27|33)l . together with the fact that fxi.P1.P2) is exponentially small (n is used 
for picking out the surviving component), we find for the normalized zero-modes at 

Z = Pm (|Cl = \j^V2/K) 

*Ux) = 'pTi^KW^ina.eUd, (p-^f . (37) 
Using the fact that (cmp. Eq. (jHJ) 

<fc\x) = 1 " 2n\C\ 2 p a e 2m<n V 3 ^, (38) 

the following linear combination recovers the bipole zero-mode 

p 2np v y a 

(39) 

By defining ^^(x) = Pa^mi{ x ) w ith p^ an orthonormal set of complex vectors, 
with p a 1 ^ = p a e 2ma ™ \( I jp we form a complete set of orthonormal zero-modes. That 
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these are solutions of Eq. (J2J) is guaranteed by the general formalism we developed, 
but one may of course check this by substitution in the Dirac equation. This requires 

didet Qdi®- 1 = 0, with $ = (l fe - 2S m t- i y\ 

For axially symmetric SU{2) solutions with charge 2 we choose p^ = £ a bPb an d 
find for the two orthonormal zero-modes at z = p 2 

* W = ^4 /2 (^) ) ^ (*), = 2*? (p«) t S- V (1) - (40) 

with 07^ = 1 — 0^ , as shown in Ref. [TT], only depending on the constituent locations 

y$ read off from the eigenvalues of Y m . Many choices of y m , £ a , a\, p a and p 2 
actually give rise to the same constituent locations, and hence the same expressions for 
Afj,(x) and ty^. It is important for consistency that this will hold for as well. Apart 
from an irrelevant phase, this is indeed the case (checked for many random choices of 
the parameters). The explicit analytic formulae in terms of the 4 arbitrary constituent 
locations, apart from the constraint on the ordering < y 2 < y^ < y 2 > read 

2 (0 _j_ W _ I (») _ (0i 

^jw = i-tf=i-n V tct <■> • (41) 

<=i r\' + r 2 > + \y\' -y 2 '\ 
07 2 Ux) = e*W^ " 2 ) + V * K 1 ^#3, 

where = |x — et/^" 7 ^ | and the constants iV, iV^ are given by 

N 



^. W-IW-M ,^ (42) 



The phase 7 vanishes when a™ = £q = 0, but is of course irrelevant for checking 
to be a properly normalized fermion zero-mode, orthogonal to iff^. In Fig. 0] we give 
an example for the behavior of these zero-mode densities. We choose y± = —6.031, 
?/2 = —2.031, y^ = 2.031 and y 2 = 6.031. These are the constituent locations also 
found in Fig. El based on the axially symmetric solution with p 2 = 1/4, a a = £0 = 0, 
£ = 3.5, Ay 2 = 1 and p% = p 2 = 2. Shown are the results for both zero-modes 
(bipole zero- mode on the right) at finite temperature, (3 = 1 (bottom), and at infinite 
temperature (top). Note that these two only differ in the cores of the constituents, 
regular at finite temperature, but singular for the self-dual Dirac monopoles one is 
left with in the high temperature limit. The bipole zero- mode density is shown on a 
scale enhanced by a factor 5. Its reduced height is due to the fact that this zero-mode 
decays much slower than the other one, as can be read off from the behavior of (j)Z\{x) 
in Eq. flUJ). 
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Crucial for the normalizability of both zero-modes is that l/0(i) is constant on the 
Dirac strings, where l/4>s vanishes 

/ (2) (2) \ 

</» (i ;(o,o,x 3 ) = i, 0^(0,0,0:3) = -^ N { 1 $Pfc ) for yr<x3<y?, 

(i ;(O,O,x 3 ) = l, 0^(O,O,x 3 )= e^iv(4^C) for yf > < x 3 < yf\ (43) 
We are now in a position to answer the question what happens when violating the 



Figure 4: The two zero- mode densities at z = //2 = 1/4 (same configuration as Fig. OJ). 
The bipole zero-mode (right) is at 5 times the vertical scale of the second zero-mode 
(left). Top for the high temperature limit and bottom for finite temperature {(3 = 1). 




Figure 5: The two zero-mode densities for Eq. (|41jl with the order of y[ and y\ 
interchanged. The second zero-mode (left) is at the same scale, whereas the bipole 
zero-mode (right) has its vertical scale magnified by a factor 10, compared to Fig. 0J 



constraint on the alternating order of oppositely charged constituents, by smoothly 
deforming from y^ < y^ to y^ > y^\ Under this deformation, both and \E f< - 2 - ) 
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remain zero-mode solutions of the Dirac equation, and remains a self-dual abelian 
gauge field. We now have a Dirac string for y^ < x% < yjj of double the usual 
strength {<p^ 1 behaving opposed to x 2 ± ), where the second zero-mode density 

diverges, as illustrated in Fig. E3 This is because 0(2) will no longer be constant on the 
double Dirac string. The bipole zero-mode, on the other hand, remains well defined. 
It actually vanishes identically on the double Dirac string (cmp. Fig. 1 of Ref. [21]), 
and no longer "sees" the two inner self-dual Dirac monopoles. 

It would be interesting if one could formulate an index theorem for these abelian field 
configurations with singularities, but this will not be straightforward as our analysis 
shows. It is yet another subtlety in describing the monopole content of non-abelian 
gauge fields. Developing a better understanding of these constraints, that affect the 
long range properties of configurations, is our main motivation for these studies. 



5 Appearance of conserved quantities 

Our analysis has shown that in all CclSGS, clS long as z stays away from the impurities, 
the zero-mode density is exponentially localized to the cores of the constituents in the 
far field limit. The sizes of the monopole cores shrink to zero in the high temperature 
limit (with masses scaling proportional with the temperature), therefore we expect 

V m (f) = (4 7 r)- 1 Tr(i? m 1 (^)), (44) 

cmp. Eqs. E]and[2ni to be harmonic almost everywhere except for singularities tracing 
the cores of type m monopoles for z G (/z m ,/i m+ i). Since the caloron gauge field does 
not depend on z, this interpretation requires V m not to depend on z. The trace in 
Eq. f)44j) is necessary to remove any z dependence due to the fact that the basis of 
zero-modes is only defined up to a - possibly z dependent - unitary transformation. 
All this is obviously true when R(z; x) is piecewize constant, as for k = 1 and for the 
class of axially symmetric solutions of Ref. [TI]. In this case the zero-mode density 
in the high temperature limit reduces to the sum of k delta function, located at the 
appropriate constituent monopole locations. 

To show that V m is independent of z, even when R(z, x) is not piecewize constant 
we solve the Riccati equation, Eq. (J2UJ) . iteratively in l/\x\ and obtain the multipole 
expansion for V m (x). We can then use the Nahm equation for Aj(z) to check if every 
moment of this expansion is independent of z. We restrict our attention to the region 
between fi m and /v+i and perform a rotation U(x) 

( sin(#) cos(<y?) sin(6 l ) sin(</?) cos(#) 



Ytiz) = Uv(x)g(z)^tgi(z), U(x) 



cos(6) cos(y) cos(6 l ) sin(ip) — sin(0) 
\ — sin(<^) cos(f) 



(45) 
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with (9,(p) defined such that U\j = Xj = Xj/\x\ (the dependence of Y(z) on x will 
always be implicitly assumed). This leaves rotations around x as a remaining freedom, 
which we will make use of later. For the axially symmetric solutions discussed in 
Sect. El Y(z) = U(x)Y m , cmp. Eq. (jMjl . 

The Nahm equation, which is invariant under rotations, is equivalent to (working 
in the gauge where A is constant, removed by the gauge transformation with g(z)) 

^-Yi(z) = —irieijk[Yj(z), Y k (z)\. (46) 

We introduced Y(z) in Eq. (|45|) such that 

R(z;x) 2 = l k \x\ 2 - 2\x\Yx(z) + Y 2 (z), (47) 
has a simple form. Writing 

R±(z) 2 = t k \x\ 2 - \x\Q±(z; If]- 1 ), (48) 
we can now formulate the Riccati equation, Eq. (|20|) . in terms of Y(z) and Q±(z; | of | ~ 1 ) , 



, , Y x (z) Y 2 (z) 1 d 



Q±(z; 


X 


- 1 ) 




X 





X 



2n\x\ dz \ 





X 


- 1 ) 




X 





tk - w r i (49) 



which can be solved by iteration, expanding in powers of 1/|5|, something that is 
easily automated. We used the algebraic program FORM [2E! for its superior memory 
management and speed to push this calculation to a high order. We find for the first 
few terms (also easily obtained by hand), 

Rt(z)/\x\ = l k -Y 1 /\x\+^Y 2 2 + Y 3 2 ±i[Y 3 ,Y 2 ])/\x\ 2 + (50) 
\{Y 2 Y\Y 2 + F 3 W 3 =F iY 2 YiY 3 ± iY 3 Y 1 Y 2 )/\x\ 3 + ... 

where the z dependence of Yj (z) is suppressed for ease of notation and any derivatives 
with respect to z are eliminated with the help of the Nahm equation, Eq. ()46|) . We 
substitute this expansion in Eq. with R m (z) = i(R^(z) +R^ n (z)), from which we 
obtain its multipole expansion. The first few terms are 

V m {x) = ^Tr (l fc + Y 1 /\x\ + i(3Y 2 - Y 2 )/\x\ 2 + lY^Y 2 - 3Y 2 )/\x\ 3 + ... 

(51) 

A number of checks can be performed on this result. First of all V m and Q± are 
invariant under any rotation among Y 2 and Y3, V m is real and Q± are hermitian, as 
should be. By construction, cmp. Eqs. (fTTl) . (|26|l and (j44j), 



j2K{x)^ a z {x) = -8 2 V m {x), ze(/i m ,// m+1 ) (52) 

a=l 
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in the far field limit, and one verifies that this integrates to k, as required by the 
proper normalization of the zero-modes. But most importantly, we verify with the 
help of the Nahm equation, Eq. (jffij), that ^V m = to the order given (whereas 
in general neither R^(z), nor Tr(i? m (z)) are constant). For the dipole term this is 
immediate, since 4-TrYi = — 27riTr[F 2 ? Y3] = 0. For the quadrupole term we have 
£Tr{3Y? -Y 2 ) = Z -±mTr{ZY x [Y 2 ,Y 3 ] - e^Y^Yu) = 0, using the cyclic property of 
the trace, etc. Note that TrAi(z) plays a special role. It corresponds to the U(l) part 
of the U (k) Nahm gauge field, and therefore decouples in the Nahm equations. As is 
obvious from the definition of R(z; x), Eq. (JHJ), it can actually be absorbed in a shift 
of x. Therefore, where this simplifies matters we may assume Y(z) to be traceless. 

Finally we check, as conjectured above, that each of the terms is harmonic. For this 
we have to note that Y(z) depends on x through the rotation U(x), see Eq. (HHJ). In the 
expression for V m (af) the x dependence is easily recovered since Y 2 is independent of 
x and Yi(z) = g(z)Aj(z)g^(z)Xj/(2T[i\x\). It is now straightforward to verify that each 
term is harmonic. When k = 1 we may use that Y x — x ■ y m is no longer a matrix, and 
one indeed finds V m (x) in Eq. (JoTj) to be the multipole expansion for (4n\x — yml)" 1 for 
that case. Since this is harmonic (for x ^ y m ), each term in the multipole expansion of 
V m (x) has to be harmonic. For arbitrary k and Y{z) we have checked these properties 
to order |a;| _14 in the multipole expansion. 

From now on we take charge 2. In this case there are 5 independent conserved 
quantities that characterize the solutions of the Nahm equation on a given interval 
z G [fi m , /i m +i], apart from the 3 translational degrees of freedom contained in TrAj(z). 
They are given by the entries of the traceless and symmetric matrix M, 

M l3 = -I (Tr(i^)i,(z)) - l5vTt(M*)M*))) ■ ( 53 ) 

One easily checks with the Nahm equation that this is conserved, as for the quadrupole 
term considered above. Although not needed here, it is known [7 that for any k a 
solution to the Nahm equation implies det(yjAj(z) — 2iriyjXj) is constant, provided 
y G C 3 , with y 2 = 0. Indeed, for k = 2 and x = 0, using that det(yjAj(z)) = 
— ±Tr(yjAj(z)) 2 , this is equivalent to M being constant if and only if y 2 = 0. Since M 
gives all the independent invariants, it should be possible to express V m (x) in terms of 
M and x. Considerable simplifications occur, because we can write 

Aj(z) = iA ja {z)g\z)r a g{z) (54) 

(absorbing the trace part in a shift of x), and use T a T\j = 5 a bt2 + i£ a bc T c to reduce the 
matrix products to scalar products, Ai a (z)Aj a (z) = M^- + \5ijA^ a {z). We do indeed 
find that V m (x) is a function of M and x only, 

V m (x) = (l + ^Mn(x) + 0(\xn) . (55) 
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where for convenience we introduced 

Mij(x) = ^ (u{x)MU-\x)).. = iTr (Y^Y^z) - \5 l3 Y\z)) . (56) 

We performed the multipole expansion for charge 2 to order |x|~ 21 . Only the odd 
orders appear, because Y(z) is assumed to be traceless. In appendix B we give the 
term of order 21, and show how from this all lower order multipole coefficients can be 
recovered. 

6 Exact results for charge 2 

— * 

In this section we will construct SU (2) charge 2 solutions for which R(z; x) is not 
piecewize constant and analyse the localization of the fermion zero-modes in the far field 
limit. We already saw in Sect. El that on each interval (/i m ,/z m+ i) we have information 
on the zero-mode density, Eq. (p)2*|) . in terms of 8 parameters. Three of these are 
associated with Try (2), which give the center of mass coordinates for the constituents 
of type m. Of the other 5, given by the 3x3 traceless symmetric matrix M, 3 are 
associated to a rotation TZ that diagonalizes M, whereas the remaining 2 parameterize 
the eigenvalues of M. They will give a scale (D) and shape (k) parameter, see below. 

6.1 Solutions to the Nahm equation 

Explicit solutions in terms of Jacobi elliptic functions were first considered in the 
context of SU{2) charge 2 monopoles j23 EB|- These solutions can be adopted for 
the calorons provided the appropriate boundary conditions, read off from Eq. (|15|). 
are implemented. In terms of the 3x3 matrix Ai a {z) defined in Eq. (}5lj) . the Nahm 
equation becomes (away from z = /i m ) 

\j- z A\z) = detiA^A-'iz), (57) 

from which we find 

~ (A{z)A\z)) = 1 3 det A{z) = lay/det (A(z)A t (z)). (58) 

The traceless part of A(z)A t (z) is therefore independent of z, once again verifying that 
M = A(z)A f (z) — il 3 Tr(^.(2;)^4*(2;)) is constant. Here we are, however, interested in 
the equation for the trace part 

-^F(z) = 4^/det (F(z)l 3 + M), F(z) = |Tr (A{z)A\z)) . (59) 
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When we diagonalize M with a suitable rotation 71, M = 7?.diag(ci, C2, c 3 )lZ t , fixing 1Z 
such that c 2 < Ci < c 3 (note that in addition c\ + c 2 + c 3 = TrM = 0), this can be cast 
in the form 

^-F{z) = ^(F(z) + c 1 )(F(z)+c 2 )(F(z)+c 3 ). (60) 
Defining F(z) + c» = \D 2 ff(Dz) and introducing D and k to parameterize the Cj, 

D 2 = i(c 3 -c 2 ), k 2 ^^^l<l, (61) 

c 3 - c 2 

shows that the solution can be written in terms of the Jacobi elliptic functions 1 

h\z) = rr, h{z) = -— , U{z) = --, k = VI - k 2 . 62 

cn k (z) cnk{z) cn k (z) 

The overall sign of the functions fi(z) is chosen such that df\(z)/dz = / 2 (^)/ 3 (^), and 
cyclic, such that in terms of these the most general solution of the Nahm equation is 
given by 

Aj(z]a,K,h,V,k) = 2mg^(z)h^(a j t 2 + VTl jb f b {Dz)T b )hg(z), V= (4tt)- 1 £), (63) 

where 1Z is the rotation that diagonalizes M and h is a global gauge rotation (leaving 
A(z)A t (z) invariant). In the Sp(l) formalism for constructing 577(2) calorons one 
requires A^{z) = Aj(—z), a property shared by fb(Dz)r b for each b. Arranging / 2 (^) 
to be odd and f\${z) to be even in z was the reason for choosing c 2 < c\ < c 3 . 

We see that k' = (i.e. k = 1) recovers the case where R(z; x) is piecewize con- 
stant, for which V^ =1 (x) = (47r|x — y|) _1 + (47r|x + y|) _1 , with y = (0,0, V) and ±y 
the two locations of the equal charge constituents (the center of mass assumed to be 
zero), see Sect. |5] The combined zero-mode density in the far field limit is given by 
—dfV m (x) = 5(x — y) + 8(x + y), the sum of two delta-functions at these constituent 
locations. Actually, two point-like constituents necessarily implies k = 1. Comparing 
3Mn/(167r 3 |a;| 3 ) = 3 Z)f=i CjX 2 /(167r 3 |x| 5 ), see Eq. with the quadrupole term for 
V^ =1 (x), (2x| — x\ — a; 2 )P 2 /(47r|a;| 5 ), one finds that c\ = c 2 , which forces k = 1. It is in 
this context that we define y = ±(0, 0, V) as would-be constituent locations even when 
1. In which way we will approach point-like constituents when k — > 1 will become 
clear in Sect. 16.31 So far we have studied the Nahm equation away from the "impuri- 
ties". We first want to convince ourselves that there are more general (not piecewize 
constant) solutions to the Nahm equation that describe calorons with non-trivial holon- 
omy, for which we need to solve the boundary conditions at the "impurities" . 

lr The Jacobi elliptic functions are defined by sn^u) = sin(ip(u)), cn^(u) — cos(</?(w)) and dn^(u) = 
v/T— k 2 sn£(it), with u = j^ u%> d0(l — k 2 sin 2 O)^ 1 / 2 . We use boldface for k to avoid confusion with 
charge k. One also encounters the notation [23 sn(u\m) for sn^u), with m = k 2 . 
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6.2 Matching at the impurities 



We will consider here the case of SU(2) charge 2 calorons in the symplectic repre- 
sentation, for which Aj(—z) = This condition is preserved under the gauge 
transformation g(z) defined in Eq. (jl3j) . but requires us to further constrain h ap- 
pearing in Eq. (jfiHj) to be generated by r 2 . It means that for z G [fix, ^2] the Nahm 
gauge field is given by Aj(z) = A, } {z; a (1) , e~^ eWT \ £> (1) , kW). For z e [/i 2 , 1 + 
we use that periodicity of Aj(z) implies Aj(l — z) = Aj(— z) = AUz), such that 
Aj(z) = Aj(z - |;a^,^,e-i e(2)T2 ,P (2) ,k( 2 )). This was studied before by Houghton 
and Kraan for trivial holonomy (p 2 — 0) |3U] , where one of the monopole types is mass- 
less. For general holonomy, the Nahm equation reduces at z — /^ 2 to (see Eq. (fT5l) ) 

A^ - I; 3 (2) , K®, e-t e(2) -, V (2 \ kM) - A> 2 ; a« e"^ (1) -, k«) = 2«pi 

( 64 ) 

At z = n\ = — H2-, using Aj(— z) = Aj(zY, the same condition is found (one can 
check [TT] that p\ = —p 2 ). 
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5 10 15 20 

Figure 6: Plotting |Ao| versus D, Eq. fU), for k = 0, 0.9, 0.99, 0.999 and 0.9999 (left 
to right). 

A particularly simple solution is obtained for p 2 = 1/4 (equal mass constituents) 
by taking the same parameters in both intervals, up to a shift a, £> (m) = V, k^ = k, 
ftM = i 3 and 6^ = 0, collapsing Eq. flMJ) to 



p 2 = Aat 2 - 2(0, 1, 0)r 2 Vk' 



(65) 



cn k (i£> 

This can be solved by taking P 2 = ±(1 2 + t 2 ), ( a = pexp(27i7a a r 2 ), with a± — a 2 = 1/4 
which gives pf = -2np 2 Pf(0, 1,0), such that 

Dk'sn^iD 



-(0,1,0) = -A3 = 7rp 2 (0,l,0). 



(66) 



27rcn k (| J D) 

For increasing D, cn^^D) will reach zero at \D = K(k), the half-period 2 . To keep 



K{k) = \l 12 d$(l - k 2 sin 2 9)- 1 / 2 , satisfying K (k) = - log(±k')(l + C(k' 2 )) [23] 
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|Aa| finite, k' has to approach zero as well. 

In Fig. El we plot |Aa| as a function of k and D, from which we confirm that, at 
fixed |Aa|, k has to approach 1 to have well separated constituents. In this point- 
like limit, monopole constituents of type 1 are located at (0, — |vrp 2 , ±V) and those of 
opposite charge (type 2) at (0, ^np 2 , ±Z>) (choosing the overall center of mass at the 
origin). This configuration has parallel magnetic moments and placing self-dual Dirac 
monopoles at the mentioned locations would give a gauge field in the bipole ansatz, 
Eq. (|35|). However, the abelian component of the gauge field coming from the above 
caloron can at best take the bipole form in the limit discussed. For k = 1 and T> finite, 
| Aa\ = 0, and one is left with two singular instantons (of zero size) at (0, 0, ±Z>). 

For this reason we now consider a class of solutions which contains regular axially 
symmetric solutions with k = 1. Again we take p 2 = 1/4, "D*" 1 ) = T>, k'-" 1 -' = k, 
( a = pexp(27rza: a T 2 ), but now with a 2 = —ol\ = vr" 1 ^, P 2 = |(1 + r 3 ), 8™ = —6^ = 9 
and 

(cos <p sin <p \ 
10. (67) 
— sin if cos (f ) 

One finds 

p 2 b = — 7rp 2 (— sin ar 3 , — sin ar 2 , T\ + cos ali) ab , (68) 
and Eq. ()64j) takes the following form 

2V(f 1 (\D) sin 8 cos if + h{\D) cos 8 sin f)r 3 + Aaii 2 

-2Vf 2 (\D)r 2 + Aa 2 l 2 
-2V{f l (\D) cosdsmf + h{\D) sin 6> cos <p)ri + Aa 3 l 2 

This can be simplified to 

Vsm{8-f){u{\D)-h{\D)) = i7rp 2 (l-sina), Ad = -np 2 cosa(0, 0, 1), (70) 
Vsm(6 + f)(f 3 (iD) + f l (iD))=iirp 2 (l + sma), Vf 2 {\D) = -\7ip 2 sma, 

It gives a three parameter family of solutions with would-be point-like constituents at 
= ((-iy'Psin(p,0,(-l) m+J Pcos(p- (-l) m ip 2 cosa). (71) 

To have an exact point-like far field limit we need to impose k = 1 , implying sin a = 
and cos^sincp = 0. The first possibility is that cos 8 = 0, for which | cos<p|P = \p 2 . 
One finds two constituents of opposite charge to coincide. Such a solution describes a 
singular (zero-size) instanton on top of a smooth caloron. Excluding this singular case 
we are left with the choice simp = 0, for which | sm9\D = \p 2 . We now find axially 
symmetric solutions with constituent locations at 

jj%> = T^p 2 ((-l) m + (-l) j | sin^r 1 ) (0, 0, 1), (72) 



= 7rp sinar 3 

= 7rp 2 sin ar 2 (69) 

= — np 2 (ri + cos al 2 ). 
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where the overall sign comes from the fact that cos a = ±1. For cos# 7^ 0, all con- 
stituents are now separated from each other, giving a regular solution. Both cases 
were already studied in Ref. [H], based on assuming that p2 is one dimensional (cmp. 
Sect. Ej). It can be shown that for 577(2) exact point-like constituents, k^ = k^ 2 ^ = 1, 
forces p2 to be one dimensional for any choice of the remaining parameters. 




Figure 7: On the left we plot k versus p for a — (k = 1), a — — 7r/100 and a = —n/2 
(giving the lower bound for k at fixed p). On the right are shown the locations, Eq. (|71jh 
of monopoles and antimonopoles (fat vs. thin curves) in the 1, 3-plane for p = 1/4, by 
varying a from — 7r (indicated by the arrows) to 0. 

Nevertheless, when sin a 7^ 0, insisting as before that equal charge constituents are 
well separated, while keeping the centers of mass of these pairs at a fixed distance 
7rp 2 cosa, one forces k — > 1 while increasing T>, and hence approximate point-like 
constituents. We will illustrate this behavior for 6 = 71/ 4. In Fig. we plot for a 
typical value of p the constituent locations as given by Eq. (fTTj) . varying a between —it 
and (given p, a and 9, one can use Eq. (fTUj) to solve for <p, V and k). We also plot k 
as a function of p for some values of <p, showing the rapid uniform approach to k = 1. 
The asymptotic behavior for a = —it/ 2 is determined by 

k' = (l + 0(k' 2 )) , D = 4V2n V (1 + 0(k' 2 )) . (73) 

It is also interesting to inspect Ai(z) in the limit k — > 1 (or T> — > 00), to understand 
to which extent we retrieve the piecewize constant behavior of R(z;x), on which the 
point-like limit is based. For this we plot fi(D(k)z) in Fig. |HJ which apart from fixed 
rotations and an overall factor D would represent the constituent locations. Since 
V(k) / (np 2 ) approaches \/2 for k — >• 1, it means we normalize the constituent locations 
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with respect to height of the jumps in the Nahm data. At the impurities (z = ±1/4) 
we therefore expect fi(D(k)z) to go to a fixed value. The plotted cases, 1 — k = 1CT 4 , 
D(k) = 15.53 (left) and 1 — k = 1CT 8 , D(k) = 33.95 (right), clearly demonstrate how 
in the bulk / 3 — > 1 and /i >2 — > 0, but that they differ at z = ±1/4, to accommodate 
the discontinuities of the Nahm equation at the impurities. The cross-over from bulk 
behavior to the impurity values scales as D~ l , and is only absent for axially symmetric 
solutions. 



0.8 - 0.8 
0.6- 0.6 
0.4 - „ 0.4 




Figure 8: Plots of fj{D{k)z) for z £ [-±, |] at k = 1 - 10" 4 (left) and k = 1 - 1(T 
(right), illustrating the approach to the point-like limit k — > 1. 



6.3 Extended structure 

When R(z; x) is not piecewize constant, it is not clear what plays the role of the 
constituent locations. We therefore do expect some extended structure for k ^ 1. For 
this reason we now come back to analysing V m {x) = (47r)~ 1 Tri?~ 1 (z) in more detail. 
Quite remarkably, the discussion in Sect. implies that the result for k = can be 
obtained from that at k = 1. For this we may use the symmetry c 2 «-> C3, which 
leaves V m (x) invariant, apart from interchanging £ 2 and £3 (M is left unchanged when 
absorbing the rotation 1Z that interchanges c 2 and C3 in U, see Eq. (J56|) ). With the 
definitions of D = AirV and k in Eq. (|6ip one finds that this implies T> 2 — > —T> 2 and 
k — > k'. Therefore, 



\f2\ \x\ 2 -V 2 + J(\x\ 2 - V 2 ) 2 + AV 2 x 2 2 

: V m (f) = ^ . (74) 

47r7(|f| 2 -D 2 ) 2 + 4P 2 x| 



For this we use that the k = 1 result, V m {x) = (47r|x — y\) 1 + (4n\x + y\) 1 , with 
y = P(0, 0, 1), can be rewritten in the form of Eq. (J74j) by changing the sign of T> 2 
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and interchanging x 2 and x$. In some sense we may say that the result for k = 
describes point-like constituents that have moved into a complex direction 3 . We note 
that Eq. (JUJ) is singular on the ring x 2 = 0, \x\ 2 = V 2 , i.e. on the circle of radius \T>\ in 
the 1, 3-plane. But there is more of a surprise: the x 2 derivative is discontinuous on the 
disk bounded by the ring (away from the disc the function is smooth and harmonic), 
V m (x) = {2ny l \Vx 2 \/(V 2 - r 2 ) 3 / 2 (l + O(o%)), with r 2 = x\ +x§. This implies indeed 
an extended structure, with singularities on the entire disk. Before showing how to 
deal with this singularity structure we consider general values of k. 

As we have seen in Sect. El TvR^-(z) is conserved as a consequence of the Nahm 
equation. Using for R^(z) the Riccati equation, Eq. (|2T)J) . the fact that Tri?~ 1 (^) = 
2Tr(i?+ (z) + _R~ (z)) _1 is independent of z imposes a severe constraint. We make use 
of this by expanding ^(2) as a Taylor series in z. The Taylor coefficients can be 
expressed in terms of the initial conditions R^(0), through the explicit solution of 
Aj(z), Eq. (|63p (we make use of the invariance under translations and rotations to 
put a = 0, h = I2 and TZ — 1 3 , as in Sect. We thus obtain the Taylor series for 
Tri?~ 1 (^), of which all coefficients should vanish except for the 0th order. This gives 
a set of algebraic equations for the initial conditions, i?^(0), encoded in and X^ 
through 

R m (0) = i(Ri(0)+R-(0)) = ViXolz+XjTj), |«(0) -R~(0)) = V(X t 2 +X jTj ), 

(75) 

Note that D enters as an overall scale factor, cmp. Eq. (jHSJ) such that 



Dependence on the coordinates and on k is mostly left implicit. The system of equations 
is of course hugely overdetermined, and some amount of good fortune was required in 
that the first 11 orders in the Taylor expansion were sufficient to solve for X^ and 
X M . Quite remarkably these imply that X — X\ — X 3 = X 2 = 0, which considerably 
simplifies the task of solving for the remaining 4 variables. We found that 

6 = X 2 -X 2 -X 2 3 -X 2 (77) 

satisfies the cubic equation 

(2 - k 2 )k 4 + 4k 2 (a; 2 - a; 2 ) - k 4 (3x 2 - x\ - x 2 ) + ((2 - k 2 )(3a; 2 - x\ + a; 2 ) - Ax 2 ) \x\ 2 
-\x | 6 - (k 4 + 2k 2 (x 2 - x\ - x 2 ) + 4x 2 + \x | 4 ) 5 + (k 2 - 2 + \x | 2 ) 5 2 + 6 s = 0, (78) 



3 This was observed before in the monopole context [27j- It is also worthwhile to point out that the 
axially symmetric monopole solutions discussed there can appear as such in the caloron context, as 
we read off for k = from Eq. I|66(l . In this case we expect to find another class of axially symmetric 
caloron solutions, with Aa giving the symmetry axes, but since we are more interested here in the 
case of well-separated monopole constituents, we did not analyse this in further detail. 
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whereas Xi/Xq, X^/Xq and X2 can be solved for in terms of 5, 



X 1 /X 
X?,/X 



\x | 4 - S 2 + 2{k' 2 (S - 2xj) + x 2 + x 2 - x 2 ) + (2 - k 2 )k 2 

4xik'k 2 

g | 4 - 5 2 + 2(6 - 2x 2 - k 'Vi - ~ ~ ( 2 ~ k ') k2 

4x 3 k 2 

|f| 4 - S 2 + 2k 2 (x 2 -x\- x\) + Ax 2 2 + k 4 



x = J_! ! v \ / ^ . (79) 

Therefore, V 2 is a rational function of 5, x and 

V2 = (5 + x 2 )(i-x 2 x - 2 -x 2 x - 2 )' (80) 

and the proper root of the cubic equation for 5 to use is fixed by V —>■ l/\x\. Indeed, 
the asymptotic expansion for S, 




reproduces the multipole expansion of V. This was verified to the 21 orders given in 
Eq. (|1(J4|) . On general grounds it can be argued, since 8 satisfies a cubic equation, that 
V 2 has to satisfy a cubic equation as well. Its coefficients (polynomials in k and x) are 
somewhat lengthy, and therefore not reproduced here, in part because we will shortly 
present an exact integral equation for V m (x), valid for any k. The exact results for 
k = and k = 1 is most easily recovered from the cubic equation for V 2 , but the 
integral representation will be valid for these two cases as well. 

Like for k = 0, we find that V m {x) is harmonic everywhere except on a disk, now 
bounded by an ellipse with major and minor axes T>, resp. k'T). On this disk the 
function vanishes, satisfying in the direction perpendicular to the disk the expansion 
V m (x) = (2n)- 1 k'\Vx 2 \/(V 2 k' 2 - x 2 -k' 2 xl) 3/2 (l + 0(xl)) (cmp. the discussion above, 
for k = 0). By introducing the "polar" coordinates (xi,Xs) = (k'r cos ip, r sin ip), in 
terms of which the ellipse at the boundary of the disk is characterized by r = V, this 
can be written as V m (x) = (27ik'y l \T>X2\/(T> 2 — r 2 ) 3 / 2 (l + 0(x\)). Care is required to 
deal correctly with the behavior at the edge of the disk when computing the Laplacian. 
We find 

_ af V m(a) = _, W ^|«^0, (82) 



with 9(T> — r) the step function, giving the following integral representation 

„ , - 1 V ft [V dr 1 
V m m = — -— H \ dip \ =o r 

2tt\x\ 4tt 2 Jo Jo^/V 2 -r 2 J( Xl _ k'r cos ip) 2 + (x 3 - r sin ip) 2 + x\ 



3) 
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Note that k', appearing in the denominator of Eq. cancels due to the change of 
variables to "polar" coordinates. By numerical evaluation, we checked this formula 
against the exact results. It gives us confidence that we interpreted the singularity 
structure correctly. Taking an arbitrary test function f(x) we find 



/(£)<9 2 V m (£)d 3 x = 2/(0) 



dip / dr . y ' — . 

Wo ^V 2 - r 2 



IT JO 



(84) 

The integral over r is well defined for any k' and can be used to check the correct 
normalization for the integrated zero- mode density, J\f(l) = 2. Eq. (}8"4*|) is also partic- 
ularly convenient for studying the limit k' — > 0. Using the fact that Af(f) is even in 
k', we may write J\f(f) = J\f (f) + k' 2 j\f 2 (f) + C(k /4 ), with in particular 



fV (f)=2f(0) 



IT J-V 



dy 



dx 



yd y f(0,0,y) 



(x 2 + y 2 )^/V 2 — x 2 — y 2 ' 



J5) 



reintroducing cartesian coordinates x = r cos if and y = rsimp. The integral over x is 
easily performed and we find for j\fo(f) 



/V (f) = 2/(0) + [ V dy sign(y)d y f(0,0,y) = /(0, 0, V) + /(0, 0, —V), 
J-v 



(86) 



whereas gets contributions from / on the line between (0,0, V) and (0,0, —V). 

Therefore, we conclude that in the limit k — > 1 two point-like constituents are found, 
and that this limit is approached in a smooth way (despite the behavior observed in 
Fig. IBJ). As an illustration we plot in Fig. El A/"(/) as a function of k 2 for a Gaussian 




Figure 9: J\f(f), Eq. (JMjl . as a function of k 2 for f(x) = exp[-10(x 2 + x\ + (x 3 - l) 2 )] 
and V = 1. 

centered at one of the would-be constituent locations, where it takes the value 1. We see 
that indeed J\f(f) reaches 1 (linear in k /2 ) for k — > 1. We also recall that Eq. (|73|) implies 
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the point-like limit is reached exponentially in the constituent monopole separation 



For any k 7^ 1 the core has an extended structure in the high temperature limit. 



the core, but suprisingly only to a disk. Note that inside this disk the singular zero- 
mode density is negative. To guarantee the proper normalization, this is compensated 
by the singular behavior at the ellipse which forms the edge of the disk. It should be 
noted, however, that the singularity structure obtained in the high temperature limit 
will only tell us to which region the core will be restricted. We will need to resolve the 
non-abelian field components inside the core to understand how this limiting behavior 
comes about. 

7 Summary and Discussion 

We studied the zero-modes of higher charge calorons, in particular in the far field 
limit. In this limit one considers the regions outside the cores of the constituents 
monopoles, where the gauge field is abelian. The constituent monopoles are most 
prominent when they all have a non-zero mass, implying a non-trivial value of the 
holonomy. Since the holonomy can be seen as a background Polyakov loop, these 
calorons are therefore more relevant for the confined phase, where the average Polyakov 
loop is non-trivial. Lattice evidence for the relevance of these configurations is steadily 
increasing ^3 EH] > and much analytic understanding of the charge 1 calorons has been 
gained j^]. These reveal n constituents for SU(n), which when well separated become 
static BPS monopoles. To localize these monopoles it is useful to consider the high 
temperature limit, for which the masses go to infinity and the non-abelian cores shrink 
to zero size. We note that the high temperature limit is equivalent to the limit of 
infinite Higgs vacuum expectation value (recall that A plays the role of the Higgs 
field in the adjoint representation). Since also chiral fermions in a caloron background 
generically have a mass, likewise going to infinity in the high temperature limit, one 
finds these zero-modes to be localized entirely to the non-abelian cores. Indeed for 
charge 1, the zero-mode density becomes a delta function at the location of one of the 
constituent monopoles. Which one of them, is uniquely determined by the holonomy, 
and the phase up to which the zero-mode is chosen to be periodic in the imaginary 
time direction. We generalize away from anti-periodic boundary conditions since this 
freedom in choosing the phase allows us to control on which constituent the zero-mode 
localizes. It thus allows us to probe the underlying gauge field configuration, see Fig. [TJ 
as has been convincingly demonstrated in Monte Carlo studies as well [TB] . 

For higher charge SU(n) calorons, k > 1, there will be kn constituent monopoles, or 



(2D). 
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more precisely k constituent monopoles of given type. Within each type the magnetic 
charge associated to the embedding in C/ n_1 (l), and the mass is fixed. In a previous 
paper jTTj we constructed axially symmetric solutions, which shared many of the prop- 
erties of the charge 1 solutions. In particular the high temperature limit gave point-like 
constituents and we have verified here that the zero-mode density indeed localizes to 
these constituents in the expected way, see Fig. 01 This, however, is not expected to 
hold for higher charge calorons in general. The reason is that one constructs these 
solutions with the help of the Nahm equation, described by a dual U(k) gauge field 
Aj(z) on the circle, with singularities at the eigenvalues /i m of the logarithm of the 
holonomy. In case of the axially symmetric solutions of Ref. JTI] this Nahm gauge field 
is actually piecewize constant (apart from some trivial U(k) gauge rotation), and its 
value on the interval z G [/i m , /i m+ i] directly determines the locations of the type m 
monopole constituents. However, in general the Nahm gauge field will not be piecewize 
constant. In that case it remained to be seen if the constituent locations are smeared 
out over all values of A(z), or worse. 

Indeed in this paper we have found for charge 2 calorons in SU(2) that in general 
the location of the constituents is smeared out over a disk. To be more precise, we 
have shown that the zero-mode density vanishes everywhere, except on a disk that 
is bounded by an ellipse. We wish to stress that this result is found by solving the 
Nahm equation on an interval, without imposing the boundary conditions. In other 
words, a disk is described by the parameters of the solution for each interval, that is 
for each monopole type, namely the center of mass, the orientation, a shape and a 
scale parameter (eight in total). These should also be the building blocks for higher 
gauge groups SU (n) since their dual descriptions differ only in the number of intervals. 
Moreover, only the boundary conditions distinguish between calorons and monopole 
solutions, and we believe our result is of value in the context of the latter as well. 

In particular we want to draw attention to the fact that (the trace of) the chiral 
fermion zero-mode density in the high temperature limit is the Laplacian of a function 
— V m (x) that is determined in terms of the Nahm gauge field and turns out to be 
a constant of the motion with respect to the Nahm equation. We have verified this 
property by computing the multipole expansion of V m (x) to a high order, but we 
expect this can be proven from the integrability of the Nahm equations, an interesting 
problem to be pursued in the future. The fact that V m (x) is conserved is a powerful 
result indeed, since it allowed us to calculate this function exactly, on which we base 
the findings mentioned above. 

This makes us conjecture that the cores of the constituent monopoles are in general 
extended, collapsing to the disk in the high temperature limit. In itself this is not 
surprising, since one knows from the study of monopoles, when two are closer together 
than the size of their core, they show an extended structure, for charge 2 indeed in the 
shape of a doughnut [2]. It should also be noted that we found, when approaching 
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the disk from above with a test function smaller than the ring, the resulting zero-mode 
density can be negative. There are two reasons why we should not be too worried about 
this. If the core would collapse to points, the proper normalization of the zero-modes 
requires the contribution from the singularity to be quantized, giving a delta function of 
unit strength. When the core is extended, we have no such constraint. In addition, our 
equations are derived by ignoring the exponentially small terms that occur in the core. 
Therefore our result is only valid outside the core, where we indeed find the zero-mode 
density to vanish. This means the cores have to lie within the disk, whereas the only 
thing we can say about their contributions is that they have to integrate to 2 (since 
— cf V m (x) adds the two zero- mode densities). Although this may seem to resemble a 
singularity like the Dirac string, it is more subtle than that, because only like-charge 
monopoles are involved here. Nevertheless, it does show that the non-abelian fields 
inside the core have an intricate behavior. 

It will therefore be interesting in the future to try and get access to the gauge field 
and to see if the field strength shows further localization within the disks we have found 
on the basis of the zero-modes. Another useful object is the Polyakov loop because it 
traces the constituent monopoles [Ej and is able to reveal extended structures in the 
context of Abelian projection jHllS21- We believe to have gained sufficient information 
to get access to the solutions at finite temperature, with a fully resolved core to answer 
these questions, if necessary by numerical means. 

However, from the physical point of view, the most important result of this paper 
is our demonstration that for well separated constituents, when the scale parameter 
D is large, the shape of the caloron zero-mode density leads to point-like constituents, 
i.e. the shape parameter k approaches 1. It does so exponentially in D and we have 
shown no structure is left on the disk bounded by the ellipse, collapsing to a line in this 
limit. Any trace left over on this line is proportional to 1 — k 2 , and therefore vanishes 
exponentially in D. This comes about through the boundary conditions A(z) has to 
satisfy at the "impurities" /i m , relating the size and shape parameters, D and k. We 
leave it to a future publication to more fully describe the moduli space of solutions, 
solving for these constraints. Also here some remarkable simplifications seem to occur, 
related to the integrability of the Nahm equations. 

The results of this paper therefore provide one further step in establishing that 
non-abelian gauge field configurations can be described on a large distance scale in 
terms of abelian monopoles. Of course it is only a small step, because we use exact 
self-dual solutions to establish these results. This has been in part because we dis- 
covered earlier ^T] , somewhat to our surprise, that it is far from trivial to write down 
superpositions of these monopole fields without having visible Dirac strings all over the 
place, that would carry too much energy for comfort. In part this is due to the crude- 
ness involved in the superposition, instead requiring a fine-tuning of the non-abelian 
tails with the exponential components in the abelian gauge field (to properly absorb 
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the return flux). It has been this problem to deal with Dirac strings that for so long 
has hampered attempts to describe an interacting theory of magnetic monopoles [3*3*] . 

In the light of this it would of course be welcome if more lattice studies are performed 
to get a handle on the dynamics of these constituent monopoles. It should be pointed 
out that instantons larger than (3 will no longer reveal themselves as lumps of size p. 
Rather there is a transition region beyond which p should be interpreted to set the inter- 
constituent distance (typically of order Tip 2 j '/?), whereas the size of the lumps is in this 
region set by the mass of the constituent monopoles. This may lead to a natural infrared 
cutoff in the size distribution of instantons, provided by the temperature (rather than 
the spatial volume). Because of this it may be worthwhile to reinvestigate the issue of 
the instanton size distribution, also with the criticism presented in Ref. [31] in mind. 



Appendix A 

We will derive in this appendix the zero-mode limit. Our starting point is the explicit 
expression for the Green's function, Eq. (JT7J). It is useful to note that this Green's 
function satisfies f x (z,z') = fl{z',z). To show that Eq. (fTTj) indeed respects this 
relation, one can use (see Eq. (|T%|) ) 

° -l k )w\zX)(X 1 ')=^(,,A (87, 

° "p^fe-nr 1 ^ ^ = (tu - F- 1 )- 1 = t 2t - (t 2t - F K )-K 

We will take z Q = p m + 0, which leads to the following decomposition of jF Mm in terms 
of contributions coming from the impurities {T m ) and from the propagation between 
the impurities (H m ), 

•?>m = T m H m -\ ■ ■ ■T2HiTig\l)H n T n H n ^i ■ ■ - T m+ iH m , 

T m = exp( 2 ® s H m = W m (p rn+1 )F m (p m+ i)W~ 1 (p m ). (88) 

Note that by definition F m (p m ) = I*.. It is convenient to absorb the algebraic contri- 
butions coming from W m in the "impurity" contributions T rn . 

(nm nm \ 
9 m + + ^j^^^W^-lW (89) 

1 i? _1 (yU ) f ^m^m) ~\~ Rm-l^m) ~t~ S m i? m (/i m ) R m _i(p m ) ~\~ S m \ 

The following identities are noteworthy 

9™± + 9™ ± = l k , 0± + — 0±_ = ±R^ n 1 ( y p m )R m _i(p m ), 2_R m (/i m )6 l ™ + = S m , (90) 
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as well as the fact that m : = W / m l 1 ( / u m )T m 1 H / m ( / u m ) can be computed explicitly 

®m = — h - J = ^m-l\^m)Rm{f J 'm) y _gm Qm J • (^1) 

Finally, introducing 

4 = ^"Vj-Tv^mQO = 8J m -! • ■ • Q 2 F 1 g\l)Q n F n _ 1 ■ ■ ■ Q m+1 F m , (92) 
we find for fi m < z' < z < /i m+ i (cmp. Eq. (JTTJ) ) 

f x (z,z') = -At: 2 ( ) V m (z)(l 2fc - ^ J^V) (°J . (93) 

It will be useful to write ^ m = F^G^+xLlfG^i-F™, with L = G w+1 F TO G m , be- 
cause K = F m _i6 m _i • • ■ Qig\l)F n Q n ■ ■ ■ Q m+2 F m+ i contains the exponential factors 
in terms of which we can take the zero-mode limit. 

With (t 2k - ^ m ) _1 = F^-Qm+xil - LKY 1 Q m+1 F m , we reduce the problem to 
approximating (l 2k — LK) -1 . For this it is convenient to write LK = LK + LK, with 

fK ++ K + _\ ~_( \ f _fL ++ 0\ f _ fO L. 



Ks {T £)• K =W- + K-h Ls U" oJ" Ls Vo l. 

(94) 

after which we find 

(l 2k - LK)' 1 = K~ x (R- 1 - L - LKK^y 1 . (95) 

As we will show next, the advantage of all this is that terms containing K±± are of the 
form K~l, K^K+-, K~ + K~l or (R—y 1 = fT__ - K-+K+\K+- and that these are 
all exponentially decaying. For the first three this is easily seen using that K is of the 
form FQFQ ■ ■ ■ FQF, whereas for the last term we recall a well-known formula for the 
inverse of a 2 x 2 matrix with as entries (k x k) matrices 

x (K ++ K + \-' f {K ++ - K + _Kz\K_ + y l (K_ + - K—K^_K + +)- x \ 
\K_+ K—J \ (K+- - K ++ Kz\K__y l (K__ - K^K^K^y 1 J ' 

(96) 

From this we find that (K ) 1 = (K x ) (hence the upper indices). With K 1 

having the form F~ 1 Q~ 1 F~ 1 Q~ 1 ■ • • F~ l Q~ l F~ l , which interchanges the role of / + and 
we conclude that K behaves as K ++ and that therefore (K ) _1 is exponentially 
decaying as well. Using 

K-i-( K £ ~ K +l K +-) KK-i-( ) (97) 
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and neglecting the exponentially decaying terms, we find the simple result 



\ \ _1 / 



or including subleading terms 

^ - LK ^ = ( -L^L-f+OiX) tTolh ) ' W 
where 0{X) = 0{K^ + K + J) + 0(K^ + K-\) + 
Using Eq. (jHEJ) and 6 m+ iF m = we find 

(100) 

This can be simplified further using L = (X x ) = (L — L h L + ^L + _) 1 (cmp. 

Eq. (pHjl l. such that 

/r(^')-(§:])V-e- +1 (2 (L i r O m i (^(^)^ i( ^- (101) 

With Eqs. (EMU), noting that Z m = -6^{6~)~ x and Z+ = (fl™)-^-* (see 
Eq. (|25)0. we find after some algebra the relatively simple result given in Eq. (|24j) . 
We note that it is not directly obvious that f* m (z,z') = f^ m (z',zY. Nevertheless, 
this is guaranteed to be true from the fact that the exact Green's function respects 
this property. All we wish to mention here, is that Eq. (JHTj) implies rather non-trivial 
relations involving f^(zY an d -R^(z)^, which could be used to explicitly verify that 



Appendix B 

Using the invariance under a one-parameter set of rotations around x for M, Eq. (|53jh 
or equivalently around (1,0,0) for M(x), Eq. (J5HJ), we can for charge 2 express the 
multipole expansion of V m {x) in the following 4 independent parameters, 

p = ^M u (x), w 2 = Ml 2 {x) + MUx), g 2 ^I(M 11 (x) + 2M 22 (a;)) 2 +^M 2 2 3(x), 

s 3 = 4M 12 (x)M 13 (x)M 23 (x) + (M x 2 2 (x) - M 2 3 (x)) (M u (x) + 2M 22 (x)) , (102) 

where |x| 2 p, |x| 4 w 2 , |x| 4 g 2 and |x| 6 s 3 can be written as monomials in x. This choice 
has the particular advantage that for charge 2 the following remarkably simple form 
can be used ^ 

V m (x) = — |^|2n+l ' 5 P a " = Ua n-U (103) 
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checked to order \x\ 21 , but likely to be true to all orders. Therefore, the result to this 
order can be read off from the / = 2n = 20 multipole coefficient 

a w = p 10 + A5p s (q 2 - 2w 2 ) + 180pV+ 315/(g 4 - 8q 2 w 2 + Aw 4 ) + 630pV(3g 2 - Aw 2 ) 
525 

+ —p 4 {2q 6 - 36q 4 w 2 + 60q 2 w 4 - 16w 6 + 3s 6 ) + 3150p 3 s 3 (g 4 - 4q 2 w 2 + 2w 4 ) 
315 

+ p 2 (5q 8 - lQ0q 6 w 2 + 560g 4 w 4 - AA8q 2 w 6 + 40g 2 s 6 + 80w 8 - 48w 2 s 6 ) 

8 

105 63 

+ ps 3 (15g 6 - 120gV + 168g 2 w 4 - 48w 6 + 2s 6 ) + — (g 10 - 50q 8 w 2 

2 8 ^ 

+300g 6 u> 4 + 5g 4 (5s 6 - 96w 6 ) + 80g 2 (3w 8 - w 2 s 6 ) + 40w 4 s 6 - 32w 10 ). (104) 
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